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SmBe has been predicted to be a Kondo Topological Insulator with topologically protected con¬ 
ducting surface states. We have studied quantitatively the electrical transport through surface states 
in high quality single crystals of SmBg. We observe a large non-local surface signal at temperatures 
lower than the bulk Kondo gap scale. Measurements and finite element simulations allow us to 
distinguish unambiguously between the contributions from different transport channels. In contrast 
to general expectations, the electrical transport properties of the surface channels was found to be 
insensitive to high magnetic fields. Local and non-local magnetoresistance measurements allowed us 
to identify definite signatures of helical spin states and strong inter-band scattering at the surface. 


Topological insulators (TI) are a new class of materi¬ 
als where gapped bulk states coexist with gapless surface 
states having linear energy-momentum dispersion rela¬ 
tion protected by time reversal symmetry [1-4]. These 
surface states are the consequence of the closing of the 
bulk gap at an interface as a result of the change of 
the non-trivial topology of the Hilbert space spanned by 
the wave function describing the insulator to a trivial 
one at the boundary of the material. The large spin- 
orbit coupling of the bulk material implies that Dirac 
fermions at the interface have helical spin polarization - 
the spin degeneracy of the Dirac fermions is lifted and 
the spin becomes transversely locked to the crystal mo¬ 
mentum [5]. This leads to interesting possibilities, chief 
among them being the existence of majorana modes in 
proximity-induced superconducting state [6]. This, and 
the obvious possible applications in spintronics has led 
the search for new robust TI materials. 

The heavy fermion SmBg has recently been predicted 
to be a topological Kondo insulator (TKI) [7, 8] . At tem¬ 
peratures below few tens of Kelvin a bulk energy gap 
opens up in this material due to hybridization of the 
conducting electrons in the 5d band and the localized 
moments in the 4f band leading to the appearance of a 
Kondo insulator state. In electrical transport measure¬ 
ments this shows up as a rapid increase in the measured 
resistance as the material is cooled down below about 
40 K. The presence of the bulk gap in SmBg down to 
atleast 2 K has also been seen through capacitance mea¬ 
surements [9] and point-contact spectroscopy [10]. Below 
4 K the resistance was seen to saturate which could not 
be explained by the standard theory of Kondo insula¬ 
tors. Recently, it has been proposed that this low tem¬ 
perature saturation of the resistance was due to the ap¬ 
pearance of topologically protected surface states in the 
Kondo gap. Existence of this surface state has been con¬ 
firmed through resistivity measurements [11-13], non¬ 
local transport [12, 13], ARPES [14-16], magneto¬ 


conductance measurements [11, 12] and point contact 
spectroscopy [10[. An important implication of the exis¬ 
tence of topologically protected surface states is the exis¬ 
tence of spin-momentum locking on the surface [5] . The 
cleanest signatures of this phenomenon comes out from 
spin-resolved ARPES measurements - till date there is no 
convincing demonstration of this for the case of SmBg. 
An alternate way of addressing this issue is through lo¬ 
cal and non-local magneto-transport studies which can 
probe the existence of spin helical surface states [17-20]. 

We have studied in detail local and non-local magneto¬ 
transport on high quality single crystal samples of SmBg. 
Our studies allow us to identify definite signatures of 
helical spin states and strong inter-band scattering at 
the surface. The samples were prepared by the floating 
zone technique using a high power xenon arc lamp image 
furnace [21] . Electrical contacts were defined on the 
(110) surface by standard electron beam lithography 
followed by Cr/Au deposition. Before lithography the 
sample surface was cleaned using concentrated HCl and 
then mirror polished to get rid of any surface contami¬ 
nants. Defining the contacts lithographically enabled us 
to determine precisely the distance between the electrical 
probes which was essential for the analysis presented 
in this letter. Electrical measurements were done using 
many different contact configurations on two different 
samples to unambiguously probe the local and non-local 
transport in the system - in the rest of the letter we 
primarily discuss the data obtained in configuration A, 
described in figure 1(a). Electrical measurements were 
performed using standard low frequency ac measurement 
techniques by current biasing the sample (10 nA-100 nA 
at T < 2 K, upto 10/r A at higher temperatures) in the 
temperature range 10 mK to 300 K and up to 16 T 
magnetic field. 

Figures 1(b) shows the measured local and non-local 
u”* voltages as a function of temperature measured at 
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Figure 1. (color online) Local and non-local transport in 
SmBg. (a) The figure shows the local and nonlocal measure¬ 
ment schemes for configureation A, arrows denote the direc¬ 
tion of current flow. Contact pads numbered 1, 4, 7 and 10 
were used as current and ground probes - these probes had 
dimensions 600 x lOO^m^. Rest of the probes had dimensions 
100 X lOO^m^ and were used as voltage probes, (b) u* (olive 
dotted line) and u"* (red solid line) measured in configura¬ 
tion A, the local voltage was measured between contact pad 5 
and 6; 3 and 2 were used for nonlocal voltage measurements. 


zero magnetic field. The local voltage v’‘ increases by five 
orders of magnitude as the device is cooled down from 
room temperature to below 1 K attesting to the high 
quality of the sample. The value of the bulk band gap Eg 
extracted from vs temperature measurements increases 
from 3.5 meV at 20 K to 5.5 meV at 9 K - this increase 
of the bulk kondo gap with decreasing temperature is 
consistent with previous reports [10, 22]. 

A significant non-local voltage appears below 15 K 
and increases rapidly by several orders of magnitude be¬ 
fore saturating at temperatures below 3 K. The temper¬ 
ature at which the non-local signal appears corresponds 
to the energy scale of the many body Kondo gap in this 
system. This appearance of the non-local signal concomi¬ 
tantly with the opening of the bulk Kondo gap strongly 
indicates that the 2-dimensional surface conduction chan- 
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Figure 2. (a) Amount of current flowing through the bulk 

states (olive dotted fine) and surface states (red solid fine) as a 
function of temperature. Below 3K almost the entire current 
flows through the surface states as the resistance of the bulk 
diverges exponentially, (b) Green dashed line and cyan solid 
line show respectively the values of Rh and Rs as function of 
temperature as extracted from our analysis, the red dotted 
line shows the experimentally measured for the device. 


nels in SmBe emerge from the bulk topological Kondo in¬ 
sulator state and are not the result of accidental surface 
states. (This was also verified by preparing the surface 
in different ways by varying the surface polishing method 
and duration of exposure to the ambient - the results were 
all quantitatively consistent with each other.) The sharp 
increase in u”* with decreasing temperature observed in 
our measurements can either be due to an increase in the 
resistance of the surface channel, Rg or an increase in 
the fraction of total current flowing through the surface 
channel. As shown later in this letter, the resistance of 
the surface channel is almost independent of temperature 
and the increase in measured u”* is because of an increase 
in surface current with decreasing temperature. 

The total current driven through the sample gets dis¬ 
tributed into three channels: 

*0=i6(T)+*i(r)+C'(T). (1) 

Here io is the total applied current bias, fb(T) is the cur¬ 
rent flowing through the bulk of the sample, is(T) is the 
local surface current and f"*(T) is the non-local surface 
current. Equation 1 combined with the fact that the 
bulk resistance is activated down to the lowest temper¬ 
ature measured [9, 10] allows us to extract the values of 
bulk and the surface current components over the entire 
temperature range - the results are plotted in figure 2(a). 
The results can be understood using the following argu¬ 
ment: At high temperatures (above 15 K) current flows 
only through the bulk channel as the bulk resistance is 











3 



0 400 800 1200 1600 2000 


x((inn) 

Figure 3. (a) Electric field profile obtained from the finite 

element simulations. The potential at contact 4 was set at 
2.94 mV and contact pad 7 was grounded. The rest of the 
pads were left floating. The arrows denote the direction of 
the tangential surface current density, (b) Plot of the line 
profile of the electric potential on the surface extracted from 
the simulations as a function of the distance x from the left 
edge of the device. The blue dashed line shows the potential 
profile on the top surface while the red dotted line shows the 
potential profile for the bottom surface. The yellow shaded 
areas mark the position of the source and drain contact pads 
on the top surface. 


very low in this temperature range. As the tempera¬ 
ture decreases, the resistance of the bulk channel Rf, di¬ 
verges exponentially due to the opening up of the bulk 
energy gap. Rg on the other hand remains almost con¬ 
stant with temperature. Consequently, the proportion 
of surface current increases rapidly with decreasing tem¬ 
perature and at very low temperatures almost the entire 
current flows through the surface channels. The exact 
division of the surface current into the local ^^(T) and 
the non-local *s^(T) components depends on the relative 
positions of the electrical contacts on the device and on 
the anisotropy of the sheet resistance tensor. Figure 2(b) 
shows a plot of the resistance of the surface states Rg and 
the bulk resistance Rb of the device extracted from our 
analysis along with the measured v^/iK The agreement 
of our calculated value with the experimental result sup¬ 
ports the picture that Rg is almost constant throughout 
the temperature range below 15 K. However, a very small 
variation of Rg with temperature is found when temper¬ 
ature in lowered below 3 K which hints at richer physics 
underlying the surface states. 

To get a quantitative understanding of the local and 
non-local surface transport in our device we performed 
extensive simulations based on finite-element method us¬ 
ing COMSOL Multiphysics modelling software [23] . The 



Figure 4. Olive filled circles show the measured for voltage 
probes separated by 200^m and located at an average distance 
d away from current probe. In this configuration current was 
sourced from contact pad 9 and pad 10 was grounded. Blue 
dotted line is obtained from the simulations for the same con¬ 
tact geometry. 


simulations were performed for the exact dimensions of 
our device. The system was modelled as an ohmic con¬ 
ducting surface with an insulating bulk. An example of 
the surface current densities and the electric field profile 
obtained for configuration A is shown in figure 3(a). (see 
supplementary section for more details on the simulation 
technique and results). The line profile of the electric 
potential on the surface was extracted from the simu¬ 
lations and is plotted in figure 3(b). The blue dashed 
line shows the potential profile on the top surface while 
the red dotted line shows the potential profile for the 
bottom surface. The yellow shaded areas mark the posi¬ 
tion of the source and drain contact pads on the top sur¬ 
face. Using this model, we calculated the dependence of 
the non-local surface potential at various distance from 
source contact d. The analysis of the non-local poten¬ 
tial was complicated by the fringing of the electrostatic 
equipotential lines as one moves away from the source 
contact. To quantify the effect of this fringing on the non¬ 
local transport, we measured between equally spaced 
voltage probes located at different distances d from the 
source probe. Figure 4 shows the measured and calcu¬ 
lated (from simulations) u"* data. The data shows that 
the non-local voltage measured between equally spaced 
voltage probes falls off rapidly as the distance between 
the voltage probes and the source probe increases. The 
excellent agreement between the measurement and sim¬ 
ulation results lends credence to our model that trans¬ 
port in the ultra-low temperature range is through sur¬ 
face states only. 

To probe the effect of magnetic field on the surface 
states in SmBg, we carried out local and non-local 
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magneto-transport measurements in magnetic fields 
upto 16 T. If the system at low temperatures indeed 
goes into a topological Kondo insulator state, then 
conventional wisdom suggests that the surface states 
should be quenched in the presence of a magnetic field. 
Indeed it has been shown in a previous work that the 
presence of magnetic impurities quenches the surface 
states in SmBg [11]. In figure 5(a) we plot the local 
magnetoresistance (MR) measured at 20 niK while the 
non-local MR is plotted in 5(b). In both cases the MR is 
very small - of the order of a few percentage even at 16 T 
magnetic field. This is in contrast to what we expect 
for TI surface states in the presence of a time reversal 
breaking field. On the other hand it has been recently 
suggested [24-27] and experimentally demonstrated [28] 
that robust topologically protected spin helical states 
can indeed exist even under the conditions of broken time 
reversal symmetry. This magnetic field insensitivity of 
topological surface states in two-dimensional topological 
insulators has been ascribed to a spin-Chern topological 
invariant. We note that the conduction through the 
surface states of time reversal invariant topological insu¬ 
lators is also expected to be insensitive to time reversal 
symmetry (TR) breaking perturbations if the Fermi 
energy of the surface Dirac cones is not at zero energy 
(which is where a gap is expected to open with TR 
breaking) and if the surface is relatively free of impurities. 

The sharp dip near zero magnetic field in the local MR, 
seen more clearly in the zoomed in graph in figure 5(c), 
is due to weak anti-localization effect (WAL) which is ex¬ 
pected in the case of systems with large spin-orbit cou¬ 
pling. The magnitude of the dip decreases with increas¬ 
ing temperature eventually vanishing above 2 K. This 
can be understood as follows: with increasing temper¬ 
ature the inelastic scattering rate of the charge carriers 
increase which eventually destroys phase coherence of the 
conducting electrons essential for WAL [29, 30] . The low 
field local magneto-condutance (j{B) at different temper¬ 
atures was fitted to the Hikami-Larkin-Nagaoka (HLN) 
equation in the presence of large spin-orbit coupling [31] : 


where is related to the phase coherence length as 
B^ = h/{4:el‘^). The fits yield a = 0.5 in the tempera¬ 
ture range 1K>T>0.1K which is the value expected 
for 2-dimensional TI having a single coherent conduct¬ 
ing channel. As the temperature decreases, the value of 
a steeply increases and saturates at temperatures below 
50 mK to a value of about 1 which is the value expected 
for two independent topological coherent channels. Al¬ 
though this doubling of the parameter a with decreasing 
temperature has been reported recently in the case of the 
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Figure 5. Magnetoresistance of the device (a) Local MR mea¬ 
sured in the range -16T to 16T. (b) non-local MR measured 
in the range -16T to 16T. In both cases the measuremets 
were done at 20mK. (c) Low field local MR showing weak 
anti-localization, (d) Low field non-local MR. 


TI material Bi 2 Te 3 [32], a proper understanding of this 
effect must await further theoretical investigations. In 
the absence of such detailed theoretical modelling we con¬ 
jecture that the saturation of the value of a to about 1 at 
the lowest temperatures indicates either the dominance 
of two surface Dirac bands from amongst the expected 
three bands or strong inter-band scattering so that there 
are effectively only two conduction channels. The value 
of obtained from the fits to the low-field MR data 
was « 700 nm at the lowest temperature. l,f, decreases 
with increasing temperature as a power law oc 
with p « 1 over the entire temperature range measured 
suggesting that the major source of scattering is electron- 
electron interactions. 

The low field behaviour of the non-local MR, shown 
in figure 5(d) is not so clear as it is in the case of local 
MR. A possible reason could be that the local magneto¬ 
transport is confined only to the top surface of the device 
and hence is decided by the dispersion relation of only the 
[110] surface. The non-local current, on the other hand 
flows through multiple (not only the top and bottom, 
but also the side) surfaces of differing crystallographic 
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orientation and hence differing dispersion relations. This 
makes the low-field non-local MR very dependent on the 
exact electric and magnetic field profiles and hence quite 
difficult to tackle analytically. Understanding the non¬ 
local magneto-transport in these systems will probably 
require significant theoretical investigations. 

To conclude, we have studied the temperature and 
magnetic field dependence of the non-local transport in 
the topological Kondo Insulator SmBg. In contrast to 
general expectations, the electrical transport properties 
of the surface channels was found to be insensitive to high 
magnetic fields. Local and non-local magnetoresistance 
measurements allowed us to identify definite signatures 
of transport through multiple Dirac bands and strong 
inter-band scattering at the surface. Understanding the 
robustness of the surface states to magnetic fields requires 
further experimental and theoretical investigations. 
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SUPPLEMENTARY MATERIAL 

CALCULATION OF SURFACE AND BULK 
CURRENT COMPONENTS 




Figure 7. Electric field profile obtained from our simulations. 


channels 


Figure 6. Different current components flowing through the 
device as a function of temperature as deduced from our anal¬ 
ysis. The values have been normalized by the bias current iq. 


At ultra-low temperatures when the bulk conductance 
has frozen out, the applied current flows entirely through 
the surface, i.e. is = io- Then, we can write 

V\T) = i\(T)Rs{T)-V'^\T) = if{T)Rs(T) (3) 

From these two equations and i\ = yig = yfo we get, 
7 =(VL)/(VL+V :;0 and i?, = {{Yl+Yt)/^o■ From the 
measured values of u"* and u* we get 7 = 0.87 , 0.77 and 
0.81 for configuration A, C and E respectively (see Ta¬ 
ble 1). Assuming that the change in the surface re¬ 
sistance Rs with temperature is negligible in compar¬ 
ison to the orders of magnitude change in i?;,, using 
= Vni(T)/Rs), we get the value of the temper¬ 
ature dependent non-local surface current. Subtracting 
this from the constant total current Iq supplied to the 
sample gives the temperature dependent local current 
i\T) flowing through the sample. Subtracting the local 
surface current f(,(T)(= Vi{T)/Rs) from the total local 
current i^{T) gives us bulk current ib{T). The values of 
the different current components obtained through this 
analysis are plotted in figure 6 as a function of tempera¬ 
ture. It can be seen that at high temperatures the entire 
current flows through the bulk of the device. At temper¬ 
atures below 3 K lb drops exponentially and the entire 
current shifts to the surface channels. The non-local cur¬ 
rent f"*, which was negligible till now, becomes a signifi¬ 
cant fraction (almost 20%) of the total current. The rest 
of the current flows through the local surface channel. 

The geometry dependent factor 7 is decided by the 
effective resistances of the local and non-local transport 


7 = i?"7(-R”'+^') (4) 

which in turn depend on the lengths and surface resis¬ 
tances of the local and non-local transport paths. R} 
depends only on the properties of the top surface as the 
local current flows only on the top surface, i?"* on the 
other hand depends on the properties of the top, bottom 
and the side surfaces. Assuming these surfaces have very 
similar specific resistances, we estimate 7 = 0.945 for our 
device. Our experiments found the 7 value to be almost 
15% lower than this value. We believe this to be due to 
two reasons. Firstly, our simulations show that there is 
a significant fringing of the electric held lines implying 
that the non-local current flows through all six sides of 
the device. The second reason is that the transport in 
this material is highly anisotropic and the assumption 
that all the surfaces have similar specific resistance is 
strictly not valid. 

To understand the experimentally obtained value of 
7 better, we performed extensive numerical simulations. 
The simulations reproduce the experimentally measured 
value of u"* and their dependence on the distance from 
source contact d as can be seen in figure 4. 

DETAILS OF SIMULATION 

The system was simulated using COMSOL Multi¬ 
physics modelling software [23]. The geometry of the 
sample was designed with the actual dimensions. The 
bulk was assumed to be an ideal insulator while the sur¬ 
face was taken to be a conducting channel of thickness 
t with resistance Rs = 2940/n. The numerical simula¬ 
tion involved the solution of the continuity equation with 
the bias applied as the boundary condition to obtain the 
potential map 7. This was done for the various conflgu- 
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Configuration 

Source 

Drain 

Local V+ 

Local V 

Non-local V+ 

Non-local V 

A 

4 

7 

5 

6 

3 

2 

C 

1 

4 

2 

3 

6 

5 

E 

7 

10 

8 

9 

6 

5 

INL 

2 

1 

variable 

variable 

variable 

variable 

2NL 

9 

10 

variable 

variable 

variable 

variable 


Table I. Description of the electrical contact configurations used in this work, the numbers refer to the elecrical probes as shown 
in figurel(a). In the case of configurations INL and 2NL, the voltage probes were varied depending on the measurement. 


rations listed in table 1 . The effective thickness of the 
conducting channel was taken to be 1 A. This is justi¬ 
fied since, it was observed that the result of simulation 
did not change by varying the effective channel thickness 
in the range t = 1-200 A(see figure 8) and the spatial 
distribution of the surface state has been reported to be 
around 10-20 A [33]. 
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Figure 8. Calculated as a function of distance of the volt¬ 
age probes from the source electrode, the numbers in the leg¬ 
end refer to the effective thickness of the conducting layer 
used for the simulations. It is seen that the various data sets 
collapse on top of each other. 























